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Abstract 

Some Mumford-Shah functionals are revisited as perturbed area functionals 
in connection with brittle damage mechanics. We find minimizers "on paper" 
for the classical Mumford-Shah functional for some particular two dimensional 
domains and boundary conditions. These solutions raise the possibility of vali- 
dating experimentally the energetic model of crack appearance. Two models of 
brittle damage and fracture are proposed after; in the one of these models the 
crack belongs to the set of integral varifolds. We have felt the necessity to start 
the paper with a preliminary section concerning classical results in equilibrium 
of a cracked elastic body reviewed in the context of Sobolev spaces with respect 
to a measure. 



1 Preliminaries 



1.1 Functions with bounded variation or deformation 

Let Q C -R" be an open set with local Lipschitzian boundary. £" or | • | denotes 
the Lebesgue measure on i?" ; Ti.^^^ is the Hausdorff n — 1-dimensional measure and 
B(r2) denotes the collection of Borelian sets from Q or the family of Borelian maps 
from n to R. For any u G L^{n, R^) let 

denote the exact representative of u ; at each point x from the Lebesgue set Qu of 
u u(x) is the approximate limit of u in x. The complementary set of Qu is Su- Su 
is a Borelian set with Lebesgue measure zero. For a given vector v G S"^^ one can 
define the approximate limit in the v direction in a point x as the number u^ix) 
which satisfies : 

lim — / I u(y) — u'^{x) \ dx = 

p^O pn J^y^Bp{x) : {y-x,v)>0} 

If a; G r^u then n^{x) exists for any u and is equal to u(a:;); denote then by Ju C Su 
the set of jump points of u , i.e. the set of points x G Su for which u'^{x) and 

vT^ix) both exists for a point-dependent v. 

The space BV(0, i?™) of i?™-valued functions of bounded variation on the open 
set O is the subset of i?"*) of all functions whose distributional derivative, 

regarded as a set measure, has bounded variation. In the same way one can define in 
the natural way the space BV^ocl-R"; R"^) of functions with locally bounded variation 
which is a subset of Lj^^{R"', R^). BV(0, i?™) is a Banach space endowed with the 
norm: 

||u||bv^ = ||u||ii+ I Du I (ft) 

where | Du \ (B) is the variation of Du over B G B(r2). 

The space BD(n) of i^^-valued functions of bounded deformation on Q is the 
subset of V-{Q,,R^) of all functions whose symmetric distributional derivative Eu 
has bounded variation. It is a Banach space too endowed with the norm 

||u||bd = ||u||li+ I EVL I {^l) 

Obviously BV(r2, i?") C BD(ri). The resemblances and the differences between 
these two spaces are visible in the following results. 

Theorem 1.1. Let u G L^i^^R""). Then 

• (De Giorgi) If u e BV(ri, i?™) then Su is countably rectifiable, 7^""^(Su \ 
Ju) = cind in l-r'~^ -almost every point x G Su exists the approximate limits 
of u in the directions i'{x) and — z^(x) where i'{x) is the normal to Su in x. 
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• (Kohn, Ambrosio, Coscia, Dal Maso) Let m = n and u G BD(ri). Let Ou be 
the Kohn set : 

eu = lxen: H^3,p l-g"l W^)) > o| 

Then 0u is countably rectifiable , Ju Q ©u o-'nd 'H"~^(©u \ Ju) = . 



Theorem 1.2. Let u G T/ien 

• ( Calderon, Zygmund) If u & BV(f2, R^) then u is approximately differ- 
entiable 

C^-a.e. in The approximate differential map x i— > Vu(x) is integrable. Du 
splits into three mutually singular measures on O 

n-l 



Du = Vu dx + [ul (8) + Cu 

where [u] is i/ie jump of u in respect with the normal direction on Su i^- Cu 
is the Cantor part of Du defined by Cu(A) = D^u{A \ Su) where D^u is the 
singular part of Du in respect to 

(Ambrosio, Coscia, Dal Maso) Let m = n and u G BD(r2). Then u has 
symmetric approximate differential e(u) CJ^-a.e. in Q and Eu splits into three 
mutually singular measures on 

Eu = e(u) dx + [u] uH?-^ + E'^u 
Moreover u is approximately differentiate C^-a.e. in Q.. 



The previous theorem allows the introduction of two new spaces. The space of 
special functions with bounded variation SBV(r2, i?™) is the subset of BV(r2, i?™) 
of all functions whose Cantor part Cu is equal to zero. The space SBD(O) of 
special functions of bounded deformation is the subset of BD(0) of all functions 
whose Cantor part E'^u is equal to zero. These special functions have the following 
regularity properties: 

Theorem 1.3 

• W^'^{Q.,R^'') C BV(r2,i?"^). The inclusion is continuous in respect with the 
Banach space topologies. If 

u G SBV(0, i?") 
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then 

u G W^^^{n\Sn,R^) 

Moreover ifu G W^''^{Q\K, R'^)nL°°{Q, R^) , where K is a closed , countably 
rectifiable set with W'^iK) < +oc, then u G SBV(J^,i?"*) and W-^{K \ 
Su) = 0. 

• Let LE^{Q) be the Banach space of L^{Q, K^) functions with L} symmetric 
differential. If u £ SBD(r2) then u G LE^{il \ Ju)- Let K be a closed , 
countably rectifiable set withW^iK) < +00. //u G LE\n\K)nL°°{n, R"") 
then u G SBD(J^) and \ Ju) = 0. 

Remind also that if u G SBD(n) then the extension of u by in i?" \ O is a 
SBDioe(i?") map. 

1.2 Sobolev spaces with respect to a measure and fractured media 

Let be a finite positive measure on fi, where Q is an open bounded subset of 
i?", with Lipschitz boundary. For any p > 1 

W^'P{n, R'') = |u G B\{n, R^): \Du\^fi , J \ dn< +oo| 

denotes the Sobolev space of order p with respect to the measure fi. Obviously 

w^f{n,R'') = w^'P{n,R'') 

Also, for any countably rectifiable surface K = K in ^ let consider the measure 
/X = + Ti^-^; then i?*^) is a subset of SBy{n,R^) and the inclusion is 

continuous.Therefore the following formula holds: 

{Du,(t))= [ Vu-(l)dx+ I [u] (8) n • dH"-i G C°°(J7, ii!"^^) (1.2.1) 

JK 

Moreover 

Wl^P{n, R^) C W^^P{^ \ K, R'') 

In [DGCL] is proved that 

is a reflexive Banach space endowed with the norm 

Du 

II llw^^.f II IIL^ ^ II ^ IIL^ 

The weak convergence, denoted by Uh u, is equivalent to: 

. Duh Du . 
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Another weak convergence, denoted by —r u, may be defined by: 

f u/, ^ u in L^(Sl,R^) 

Let consider the decomposition /x = aCP" + /x* and the functional 

F{x) = [ /(— )a dx 
where / is a quasiconvex function, i.e. 

fiz)< [ fiz + D(j){y)) dy , Vz G ii:"^'^ and V<^ G i?*^) (1.2.2) 

J[o,i]" 

The following result has been proved in [AFB] : 

Theorem 2.1. ( Ambrosio, Buttazzo, Fonseca) -.Let p > 1 and f a quasiconvex 
function such that 

< fiz) < C{1+ I z \P) \fz e 

If ^ € then the functional F is seq. lower semicontinuous on W^'P{0,, R^) with 
respect to the convergence —r. 

If / is a convex function then F is convex too. In the case /x = £" + if 

W'-^{K) < +00 then for any sequence {uh)h C W^'P{n,R'^) bounded in the L°° 
norm the convergence implies the convergence — Therefore F is seq. lower 
semicontinuous on every set 

Wl,^P{n,R^,M) = |u G Wl,^P{n,R^) : ||u||oo < m} 

This is sufficient to prove the existence of a minimizer of F over the subspace 

12 L 

Wn' {Q,R ,uo) of all u with the prescribed value Uq on dO. ( uq in L°°). The 
coercivity of F is equivalent to: 

lim -P'(u) = +00 

In the simplest case, corresponding to the Dirichlet problem for a fractured elastic 
body with natural configuration Q,\K , n = k, f has the form: 

f{z) = \cz:z 

and the property 

f{z) > c I z''"" f \/z G i?"^" 
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where C is the hnear elasticity 4-order tensor, i.e. it has the properties: 

^ijkl ~ ^jikl ~ ^ijlk ~ ^klij 

In this case (and with the hypotheses on fi and K) the following inequality holds: 

Vu G W^^'P^n, i^^ M) ||u||vi.i,2 < ||u||^;^,M < ||u||vi.i,2 + C{M, K) 

So F is coercive whenever the same F defined over W^'''^{Q. \ K, R!^) is coercive. The 
conclusion is: 

Theorem 2.2. : Let ft be an open bounded subset of RP' , with Lipschitz boundary 

and K = K d a countably rectifiable hypersurface with Tl^~^{K) < +00. Let 
li = C^ + n^-^ and Uo G W^'^(J), R") n i?"). Then the functional: 

F:W^^^{n,R\uo)^R 

F{u) = \ f CVu : Vu dx , 

where C is a symmetric positive definite 4-order elasticity tensor, has a minimizer, 
unique to an element of 

Rl^\n,R'\0) = {u G W;,'\n,R^O) : e(u) = o} 

Any minimizer u has the properties: 

• u G Wff^ifl \ K, R") Vp > 2 

• divCVu = in Q in the sense of distributions 

• (CVu) n = on K 

• V(/)G C°°(J),i?"^") 

/ divcp -u dx+ j Vu:<j) dx+ f [u] • {<j)Ti) dK"'^ = [ (cpn) ■ uq dH 
Jn Jn JK Jon 
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The uniqueness of the minimizer u follows from the expression of F in a classical 
way. The properties of u come from the Euler equation: 

/ CVu : Vv dx = Vv G VF^'^(0, i?", 0) 
Jn 

and from the remarks from the beginning of the section. 

The kinematic of a fractured elastic medium is therefore well settled in the 
frame of the Sobolev space W/,'\fl,R''). The minimizer u is a physical solution of 
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the Dirichlet problem since the formula ()1.2.1() shows that there are no forces or 
moments concentrated at the edge of the crack K (the border of K). 
The space of admissible stresses is (according to [DP]) : 

Yk{^) = {ct = ct"^ G L2(17, : diva = , cm = on 

Here an means the trace of o" on multiplied by the normal n. The definition is 
good because of the following proposition: 

Proposition 2.3. (Del Piero) : Let Y{Q) he the space: 

Y{^) = {cj = cj^ G L2(17,i?"X") ; diva G L'^{Q,R'') , } 

and K any smooth surface in Vt. Denote by [an] the jump of an over K . Then 
[an] = 0. 

The machinery of convex duality (see [M]) holds in this case. Let consider the 
following separate bilinear form: 

{a,e) = / a : e dx 
Jn 

The polar of the convex functional 

W : L\n,Rl^^ ^ R W{e) = \j^Ce:edx 

is defined by 

W*:L\n,R^^^^^R^{+^} 
W*{a)=su^{{a,e)-W{e): e e L\n,R^,^:i} = \ [ C'^a : a dx 
The following inequality holds: 

W*{a') + W{e") > (a',e") 

and the equality is attained whenever a' = Ce" . If u denotes a minimizer of F, then 
CVu G Yk{^)- 

For any a G Yk(17) and for any v G Wl^'^^Q, ii") n W^'P{n \ K, i?") 

(a, e(v))= / aijVijdx+ / [ui]{an)i dH"^'^ = {Du, a) = 
Jn J K 

= - [ Ui{diva)i dx+ [ (an) ■ uq dT^""^ = [ (an) ■ uq dTi"^^ 
Jn Jan Jan 

This, together with the inequality involving W and W* , give the following 
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Proposition 2.4. : For any a € Yk{^) and for any minimizer of F the following 
inequality holds: 

[ (an) • uo rfH"-^ - J / C" V -.adxK F{u) 
Jan 2 7q 



2 Introduction 

Any Mumford-Shah functional is related to a perturbed area functional. The 
result of the first section shows that, in particular cases which model the fiber-matrix 
decohesion in some composite materials, the discontinuity set of the minimizer of 
the classical Mumford-Shah functional minimizes the area functional too. 

The decohesion between the fibers and the matrix in a composite material can be 
modeled as fracture of the matrix. In the case of a composite with rigid fibers and 
elastic matrix the fracture is brittle and it has the peculiar property that a crack 
may appear on the boundary of the matrix. If this happens then the matrix no 
longer satisfies the boundary displacement imposed by the fibers and the decohesion 
appears. The problem is therefore to model the appearance of a crack in an clastic 
body submitted to a given boundary displacement. The purpose of the first section 
is to show that in the case of antiplane displacement the energetic approach to the 
problem can model the decohesion, i.e. the crack predicted by the model may lie 
on the boundary of the body. We prove that in a particular class of geometries 
and boundary conditions the only crack that may appear is a geodesic (i.e. length 
minimizing set). The result is quantitative and experiments can be made in order 
to validate the energetic model. 

In the second section we consider two phenomena that occur in the damage of 
an elastic body: brittle (or quasi-brittle) damage and brittle fracture. These two 
damage mechanisms are the macroscopic effects of the appearance of a crack. That 
crack may be formed by a macro-crack and a large number of meso-cracks; the 
existence of the latter explains the weakening of the clastic macro-properties of the 
body. Brittle damage is modeled as something occurring in the volume of the body 
and the theory predicts the existence of fine mixtures between the damaged and 
sound states (see [FMu] and [FMa]). This is a mathematical consequence of the fact 
that for a minimizing sequence {uh, A^) of the total energy functional 

E{v,A)= / w(yu,A)dx + -f / XAdx (2.0.1) 
Jn Jn 

(where A is the damaged region) the sequence (XA^) might not converge to a char- 
acteristic function of a set, but to an / : i7 ^ [0, 1]. The reason is that there is no 
control of the perimeter of the damage region. 

The energetic model of crack appearance (in mode 3, antiplane) predicts a pri- 
ori non-smooth fractures. There are three fundamental facts about the cracks K 
predicted by this model: 
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• 1. K is not a curve that fills a volume: ?Y" ^[K \ K) = 

• 2. there is a piecewise smooth manifold containing K\S where S \s a closed 
H"-i-negUgible set (see [DS], [AP] and [AFP]). 

• 3. the curvature of K is proportional to the jump of the bulk energy density 
[w\ (this relation is true only in a formal sense due to the lack of smoothness 
of K) 

Our purpose is to couple brittle (brutal and total) damage with brittle fracture. 
This can be done in two ways; the first is to add to the crack K the essential frontier 
of the damaged zone and add the volume of the damaged zone in the expression of 
the energy functional. The surface term in the expression of the energy will control 
then the perimeter of the damaged zone as well as the length of the crack. In this 
model no fine mixtures appear. 

The second way to couple brittle damage with brittle fracture is to renounce 
to control the perimeter of the damage zone; then, in the portions of the crack 
embedded in sound regions (because now fine mixtures may appear), the curvature 
of the crack has an upper bound. This is due to the fact that the volume of the 
damaged zone controls the jump of the bulk density hence the curvature of the crack. 
In this model we choose to extend the class of n — 1 smooth manifolds to the class 
of integral varifolds. 

3 Decohesion in composites and Mumford-Shah func- 
tionals 

The idea of the energetic models of crack appearance is to enlarge the class of 
admissible displacements of the body from the space W^''^{yt,R^) to SBV(J7, i?") 
or even to SBD(r2, i?"). The jump of the displacement (let say from SBV(r2, -R")) 
represents a crack; in this way all the displacements corresponding to some arbi- 
trary crack in the body lie in SBV(ri,i?"). In the preliminaries we see that the 
correct setting of the equilibrium Dirichlet problem for an clastic cracked body is in 
WI''^{VL, i?"), with 11 = U'+rC^-^ {K is the crack); this is a subspace of SBV(0, R"). 
Suppose that on a non-negligible part r„ of the boundary of the body 

dn = t;^uTJ ,Tu nr f = ID ,0 <n''-\Tu)-n''~\T f) <+oo 

the boundary displacement uq is imposed and is stress free. For technical reasons 

uoGSBVioc(i?",i?")nL°°(i?",i?") , H"-\F/\Suo) =0 

Here Su denotes the jump set of u; because u G SBV/oc(^"', R"") the jump set 
which is usually denoted by Ju coincides modulo a Ti"^ ^-negligible set with the 
complement of the Lebesgue set Su- The space of admissible displacements is 

SB\{n, i?", uo) = {u G SBV/oc(i?", i?") : u = Uq in i?" \ H} 
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We are searching for a minimizer u G SBV(0, i?", uq) of the functional / over the 
same space, where the total energy 7 is a Mumford-Shah type functional: 

/(v) = 1 [ CVv : Vv dx + G'H''~\S^^) 
2 Jn 

We remark that / is well defined over SBD(r2) if one replaces Su with Ju- The 
existence of a minimizer of I in this space is proved in [BCDM]. 

In the case of antiplane displacements — n = 3 and u(,ti, ,T2, 2:3) = (0, 0, u(xi, X2)) 
— the correct setting of the problem is to consider O, G and u G SBV(ri, R, Uq). 
The expression of the total energy / becomes 

7(v) = I Vv p dx + GH'^-i(Sv) 

It is known from [Amb] and [DGCL] that the latter / is Lj^^ lower compact so 
there is at least a minimizer u of / in the space SBV(0, i?, uq). Much is known 
about the geometric properties of the crack predicted by this model, i.e. Su; among 
them we mention the fact that Su is a subset of a smooth hypersurface, excepting a 
TC^- negligible set (recently proved in [DS]). That push us to return to the original 
Mumford-Shah functional and minimization problem. Consider the functional 

J{v,K) = 11 CVv : Vv dx + GTiJ'-^{K \Tf) 

defined over the set M(uo) of all (v, K) such that K = K <Z Q \s a piecewise 
hypersurface, C K and v G -R"), where /i = £" + H"^^^ . Also v = uq on 

r„ \ K. With the help of J(-, •) we can define a perturbed area functional over all 
piecewise closed hyper surf aces: 

J{K) = m.ia{J{w,K) : {^^,K) G M) 

J(-) may or may not have a minimizer and this is the reason for using the functional 
I first proposed by Ambrosio. However it seems natural to us that for a big G the 
minimizer (if it exists) of J(-) minimizes also the area 7^"~^(-). Such a situation is 
described in the sequel. 



For any smooth K (i.e. with the properties listed several times before) denote by 
uk the displacement of the fractured body Q.\K under the imposed displacement 

Uq. 

J{VLK,K) = J{K) 

We restrict to the case of antiplane displacements so u is a scalar function and 
the bulk energy is a quadratic expression in Vu. Suppose that ur^ G and 
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< c <| Vupy. |< C everywhere in il. Define then u to be any function with the 
properties: 

du _ 1 dvLTf 
dxi ||Uo||oo 9X2 

du _ 1 duTf 

9x2 lluolloo 9xi 

The level surfaces of u form a congruence of curves in The free surface is 
tangent to the congruence. This congruence defines a system of open (smooth) 
neighbourhoods V{u) by: 

\/A € V{u) dA \ Tu is locally a level set of u 

For any A € ^(u) we will denote by 9uA the part of the boundary of A made 
by curves u = ct. Let be a piecewise smooth curve and O' G ^(u) such that 
K\rf cn'U du^'- Recall that 



J{K) = J / I VuK f dx + Gn\K) 



The following field is an admissible stress (i.e. it belongs to Yxiurf), see the pre- 
liminaries) : 

Vur^ n\n' 
n' 



Proposition 3.4 (which affirms the familiar minimum principle in stress) gives 
J / I VuK P dx > / (an) • uq dH^ -]- f \ a \^ dx 

^(r/) - J{K) <-\f\ Vur, |2 + / _(Vur,n) • uq - Gn\K \ T f) 
^ Jn' JTunn' 

Denote by r the clockwise tangent to dCl'. After some work we find 

J(^f)-J{K)<h\no\\oo f ur, ■{VuT)-Gn\K\Tf) 

Suppose now that r„ has two connected components T^ and and uq is taken like 
this: 

Uo = on and Uq = A = ||uo||oo on 
Consider now the projection, called P, along the congruence u = ct. on F^. Then 

1 / /Vur.n) .udn'> I A' Vu{P{K)) 

2 Jd^n' 2 
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where Vu{M) is the variation of u over the set M, n\M) < +00. Moreover 
inf i J / (Vur.n) ■ u dH^ : n' e V{u), K \ T/ C fi' I = ^A^ Vu{P{K)) 

[ ^ Jdun' J ^ 

Remark also that 

Vu{K\Tf) > Vu{P{K)) 

Take D{K) to be the intersection of all members of V{u) containing K \ T f and 
denote by T{K) any of the smallest geodesies in D{K) which separate the two 
connected components of dD{K) D T^. Then 

n^{T{K)) < n\P{K)) , 7i\T{K)) < n\K\rf) , Vu{P{K)) < Vu{T{K)) 

The following inequality holds because on the assumption on | Vur^ |: 

^AVu(K) - Gn\K \ Tf) < H\K \Tf) Qa^C - 

Therefore if 
then 

JiTf) - J{K) < \a^Vu{K) - Gn\K \ r^) < 
Also, remark that 

supi^Vu{P{K))-GH^{K\Tf) : K = K^ = snv^Vu{T{K)) - Gn^{T{K)) : K = K^ 
\/\'Vn{V{K)) - Gn\T{K)) > Vu{r{K)) Qa^ - G^) 

Therefore if 

A^ > — 

c 

then the following measure over the a-algebra over ^(u) is a positive one: 
n' ^- I \ Vur^ p dx - GH^{T{^^x'D:)) 

and 

J{Tf) - J{K) <l [ I Vur, P dx - Gn\r{duD{K))) 

2 J DiK) 

The maximum of the right term from the previous inequality is attained in D{K) = 
ri and for any smallest geodesic in O U dyi^ separating the two components of T^ 
this inequality becomes an equality. 

Proposition 1 : With the assumptions on Vt and Uq mentioned before, there 
exist two positive constants m < M , depending only on O, T^, G, such that: 
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• if 



|uo||L < rn 



then J(rf) < J{K) for all admissible K; this means that no crack appears 
under the imposed antiplane displacement Uq, 

if 

then any smallest geodesic in O, which separates from is a minimizer of 
J. 



The constants m and M can be precisely determined in some particular cases. 
For example if = 5(0, R)\B{Q, r) is a ring and Uq = on dB{0, R) and Uq = const. 
on dB{Q, r) then 

m = M = r In — 
r 

The only crack that may appear in this case is dB(0,r). 

If O = (0, a) X (0, L) is a rectangle with the free surfaces {0} x (0, L), {a} x (0, L) 
and Uo = on (0, a) x {0}, Uq = const, on (O.a) x {L}, then (calculus done by J.-J. 
Marigo) 

m = M = L 

The only crack that may appear in this case is any line (0, a) x {.x}, with x € [0, L]. 

There are two types of mechanical experiments that can be made in these cases 
in order to validate the model. One may try first to compare, for a given brittle ma- 
terial and geometry, the Griffith material constant with the constant G determined 
from the critical displacement measurements. Second, for the same material but for 
different geometries, the experimental values of G can be compared in order to see 
if G is a material constant. 



4 Perturbed area functionals in brittle damage mechan- 
ics 

Since the Mumford-Shah functional seems to be intimately connected with the 
area functional then the methods and objects used in the variational problems in- 
volving the area functional may be useful in some free boundary problems. We 
present in this section two different models of brittle damage coupled with brittle 
fracture. In thc^ first we use the SBD approach but in the second we enlarge the 
class of n — 1 smooth manifolds to the class of integral varifolds. This is not a 
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new idea; in the [AFP] and [AP] the authors use varifold techniques to prove the 
regularity of the minimizer of the Mumford-Shah functional. 



Let consider the following space: 

M = {{u,A) : V{A,n) < +00 , u G SBD(n) , E (uxa-) < Exa} 

We define on A4 the energy functional {A is the sound region; from here the minus 
sign in the expression of the energy) : 

J(u, A)= [ [w{eiu)) - 7] dx + GT-r~^ (Ju U d* A) 

J A 

Theorem 1 Let (u/^, Ah) C Ai be a sequence such that for a// /i € N 

{\\uh\\oc + Jiuh,Ah)}<C (4.0.1) 

Then one can find a subsequence (u/i^,^/i^) and an element (u, ^) € M such that: 

u/,^ u strong in ij^^i^, ii") 

XA,^ ^XA in BY{n,R) 
J{u,A) < lim inf J (u/,^ , A^^ ) 

Before giving the proof, remark that this theorem gives the existence of the 
minimizers of J in the set 

M{uo) = {{vL,A) : (u|„,A) eM ,uxA = no on 

where uq € SBD;oc(-R") H L°°{R'^,R'^). Remark also that for any minimizing pair 
(u,^) ||uxyi||oo < ||uo||oo- By a translation argument one can prove then that the 
functional 

K{u,A)= [ [w{e{uxA)) - IXa] dx + GW-HJuxa) 
has minimizers in 7W(uo). 

Proof of the Theorem 1. By the compactness theorem for sets with finite perime- 
ter it follows that we can extract a subsequence of (u^, A^) (denoted by the same 
name as the initial sequence) such that 

XAh ^ XA 

in BV(r2,i?). Since {uh,Ah) € A4 then u is a linear antisymmetric function out of 
Afi. The assumptions on w, (|4.0.1j) and £"($7) < +00 imply that for all h: 

[ \uh\ dx+\ E^Uh im+c [ I eiuh) |2 dx + 7^"-HJuJ < C 
Jn Jn 
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By the compactness theorem proved m [BCDM] we can extract a subsequence such 
that (keeping the notation): 

u/j ^ u strong in Lj^^{n, i?") 

E^Ufi E^u weak * as measures 
eK) ^ e(u) weak in ^1(17, fi^^^^^) 
Therefore Euh Eu weak * as measures and XA^ XA" in eBV{0,, R), hence 

E (^UhXAi) E {uxA'=) weak * as measures 

()4.U.1() imphes again that the densities ^Exa^''^ uniformly bounded in L°°, 
therefore one of the Hutchinson theorems (see [Hut], 4.4.2.) assures us that 

E{uxA<^) < ExA 

The lower semicontinuity of J is a straightforward consequence of the previous facts. 

Recall that a (n — 1) varifold is a positive Radon measure on the bundle x 
G{n, n — 1) with the grassmanian 

G{n, n — 1) = {S = (Sjj) : Sij = 5ij — rijiij , n € 5"} 

as fiber in any x G 17. The weight of the varifold V is 

\\V\\{A) = V ({(x, S) : X e A}) \/A e B{n) 

V^^^(S) denotes the density of V in the fiber of x € O of the mentioned bundle. 
The class of unoriented n — 1 integer varifolds is denoted by IV {i}, n — 1). One can 
consider also the class of oriented n — 1 varifolds (Radon measures over the bundle 
X G°{n, n — 1)); if q is the standard map from G°{n, n — 1) to G(n, n — 1) which 
associates to any oriented hyperplane the unoriented one, then maps oriented 
varifolds onto unoriented ones. The class of integer oriented n — 1 varifolds is denoted 
by IV°{n,n- 1). 

The first variation of the varifold V is denoted by (5V (with the total variation 
as measure | (5V |) and it has the definition: for any g G Cg(ii", i?") 



s^ig) = J giA^)Sij dv(x,5) 



In the following (u. A, K) is an admissible triple if u G W/,'\n,K^), /i = + 
n]'~\K CA is a piecewise smooth surface and XA G BV(r2, R) open with smooth 
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frontier. A denotes the sound region of the body 0, (from here the minus sign in the 
further expression of the total energy). Any minimizer (if exists) of the functional I 

I{u,A,K)= [ {w{Vn)-^)xA dx + G7e-\K) 
Jn 

in the set of all admissible triples with u satisfying a Dirichlet boundary condition is 
a critical admissible triple in the sense that for any tj G W^''P{n, fi"), supp rj compact 
in n, 

*,/(u, A A-) = lim inf '(-■4'7'-MAU(K))-H..A,K) ^ „ 

' e^o I e I ~ 

where is the one-parameter group generated by rj. Here w is a strict convex 
positive function with the following properties: 

• w{B) = wiB'y^) for all B e i?"^", 

• there are two positive constants c and C such that c [ B^^"^ p< wi^B) < C \ 

Qsym |2 j^i^ Q g ^nxn_ 

The first variation of / with respect to rj is 

6,jl{u, A,K) = / (?i;(Vu) - 7) divrj - t^Vu Vr/ XA dx + G diVsV dH"~^ 

Jni dV J JKudA 

(4.0.2) 

where divs is the tangent divergence with respect to a surface with normal n: 

diVsV = Vhi - ViJ^i^j 
In fact if we see K as a varifold then 



SK{r]) = [ diVsV dW"~^ 
Jk 



Denote hy H = —diVsV {v is the normal to K) the scalar mean curvature of 
K. If (u, A, K) is a critical point of / and if (u, A, K) is balanced (i.e. o"(u) = 
w v(Vu) G Yk{A)) then 

[w(Vu)] = GH 

on K in the distribution sense. Also, a simple integration by parts in (|4.0.2|) gives 
for balanced triples: 

5r,I{vi,A,K)= [ div[{w{Vu) -j)r]] dx (4.0.3) 
Jn 

For an uq S SBV,oc(-R", i?") n L°°(i?", i?") a triple (u, A, K) is UQ-admissible if 
it satisfies the supplementary condition: 

u G SBV(rj,i?",uo) 
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We shall work with the set T(uo) of all UQ-admissible triples (u, A, K) satisfying 
the following property: 

3iV = iV(u, A,K)c,R Mie C°(0, [0, +00)) N f i ATC'^ + / diu^^n dH"~^ > 

JK Jk 

For any element {u,A,K) G T{uq) K has bounded scalar mean curvature, but this 
bound depends on the element. For any positive number N we define the space 
T(uo, N) to be the set of all (u, ^4, K) € T{\iq) such that: 

€ C°(l^, [0,+oo)) N j idn''-^+ [ divs^n dn''-^ > 
Jk Jk 

Off course 

r(uo)= U r(uo,A^) 

For any (u, A, K) € T'(uo, A^), if there is a subset D oiQ. such that /j^^^ w{V\i) — 
7 dx > then 

I{\i,A,K) > I{u,AuD,K) 
so the minimizing sequences have to be searched in the set of admissible triples with 

w(Vu) - 7 dx < (4.0.4) 



Ib, 



'Br(x) 

for jC^-all X G and for all r > small enough. 

If we keep {A, K) and vary u we always find a minimizer of / in this class. This 
minimizer is obviously balanced. Therefore we have two simple operations that lower 
the value of I. The question is: for a given and a given {m,A,K) G T(uo,A^), 
is there a balanced triple in the closure of a r(uo,M) with respect to a particular 
convergence, satisfying H4.U.4|) . with lower energy I? 

Suppose that we repeat the "cut" and "balance" operations without finding the 
desired triple after a finite number of iterations. We obtain a sequence {Mh^Ah,K). 
For any h G ||u/i||oo < ||uo||oo- Therefore, up to a subsequence, u^XA^ converges 
weak * in L°° to a u. 

Consider now the sequence of measure-functions (see [Hut] for definition and 
properties of measure-functions) (xa^-C", e(u/i))ft,. It is easy to see that there is a 
constant M with 

/ w{e{\ih) XAy, dx < M 
Jq. 

for all h G N . Then the compactness theorem 4.4.2. [Hut] implies that, up to a 
subsequence, 

e{vih)XAh : (/) dx ^ / e : 4> dfi 
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for all (j) G C^{Q, R^y^), where (/U,e) is a measure-function and fi <C G [0, 1] 

a.e. in Also, 



' u;(e) dfi < liminf / w{e{uh) dx 



Suppose that w is the well-known C°° energy potential of a linear hyperelastic 
material with the Hooke tensor C. The weak convergence of the sequence of measure- 
functions implies that for all v G C^(r2 \ iT, i?") 

/ Ce(v) : e{uh)xAh dx Ce(v) : e 
Jn Jn 

Recall that any u/j is balanced and C is a symmetric 4-order tensor; therefore for 
all V G C^{n\K,R'') 

/ Ce : e(v) dfi = 
Jn 

which means that 

■^Ce G YK{n) (4.0.5) 
For any h G N and for any G C^{^, R^ym) the following equality holds: 

/ Uhdiv(j)XAh dx+ / (j): e{uh)xAh dx + [uh] n : dW~^ = 

Jn Jn JxnAh 

Up to a subsequence ([u/j] nxAh)h converges weak * in L°° to a and the limit of 
the previous inequality is: for all (p G C^(r2,-R"^") 

/ udivcj) dx+ / e:</)d/i+ / a : (j) dW~^ = (4.0.6) 
Jn Jn Jk 

That justifies the notation e = ^^{vl). ()4. 0.5(1 shows nothing but the fact that u is 
balanced. 

This fact suggests to enlarge the domain of admissible triples in order to contain 
triples like {u.,^^K) which satisfy H4.U.6|) . /i <C G [0, 1]. The enlarged energy 

functional is: 

I{vi,^x,K)= [ w{e^{vi))-^ d^ji + G7e-\K) 
Jn 

We shall enlarge even more the set of admissible triples by allowing K to be an 
integral varifold. With the same kind of argumentation as before we can prove that 
there exist a minimum of the functional 



/(u,/x,V)= / «;(e^(u))-7d^ + G||V| 
Jn 



in the set of all admissible triples such that the first variation | 5V | is uniformly 
bounded. We do not repeat the proof because the use of the compactness theorem 
for integral varifolds is the only new element that appears. 
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Denote by X{\iq,N) the closure of r(uo, A^) with respect to the convergence: 

{uh,Hh,^h) (u,/i, V) iff 

—5- (u, /x) as measure-functions 
V/i ^ V as integral varifolds 
We have proved that for any sequence (u^,/x^, V^) G X{uq,N) such that 

is uniformly bounded, there is a subsequence converging to an element (u, /i, V) of 
a ^(uq, M) such that 

I(u, fi, V) < lim inf I{uh, Hh, ^h) 
Also, for any (u, A, K) G T(uo, iV) exists M e R and (v, /x, /sT) e X(uo, A^) with 

||v||oo +-f(v, /U, V) < +00 

and V balanced, i.e. for all w E C^{U\K, K"-) 

I Ce^(v) : e(w) d^ = , 
Jn 

such that 

lk(e;.(v))-^i|oo<7 (4.0.7) 
and /(u, A, K) > /(v, ^u, X). In order to prove the existence of a minimizer of / over 

X(uo) = U X(uo,iV) 

it is therefore sufficient to consider minimizing sequences formed by balanced terms. 
We use (|4.0.;jp to show that for 

any minimizer of I over X(uo, M) belongs to X(uo, A^). 

Consider a sequence of balanced {yih, A^^Kh) G ^(uo) converging in X(uo) to 
(u, /X, V). Suppose more that H4.U.7|) is true for all {u.h.,Ah.,Kh) and I{uh, AhKh) is 
uniformly bounded. Since 6r^I{u, V) can not be defined for any element of X{uq), 
we adopt the following definition of the first variation of /: take a ?? G Cq{Q, R^), 
I r] \< 1, and denote by (j)^ the one-group parameter generated by rj. For any 
sequence {uh, Ah, K^) in T{uq) converging to (u, ^,V) it is not hard to see that 
{uh.(l)~^,(p^{Ah),(j)^{Kh)) converges up to a subsequence to an element of X{uq) 
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denoted by ^,(j)^{fj,),(j)^{\')) (attention, this element depends on the sequence 

{uh,Ah,Kh)). By definition 

5,Iiu, V) = mf {hminf ^("■^^'^ |-(^)) " ^) : (u., A, /^.) - (u, ^, V)} 

(|4.U.3() and H4.U.7|) show that (5^/(u/i, Ah, Kh) is uniformly bounded. For any (u, A, K) G 
r(uo) the following function is lower semicontinuous: 

5„I{\i,A,K){e) = { 9\ ^ 

^ \ V(u,yl,i^) if e = 

Again the compactness theorem for measure-functions applied to the sequence 

[u;(6(u,)]n,) 

associated to (u/^, A^, K^) affirms that (up to a subsequence) {7iV'~^ , [w(e(u/i)]n/i) — > 

(||V||,A) with ||A||oo < 7- Because of the uniform upper bound of 6ril{uh, A^, K^), 
the mentioned lower semicontinuity implies that 

6^1{u,n,\')< [ 7? • A d|| V|| + G(5V(7?) 
Jn 

It is obvious that for any minimizer of / over a X(uo,M) 

V(u,/i,V) >0 
The latter two inequalities end the proof. 

Theorem 2: The functional I has a minimizer over the set X{uq). Moreover, 
any minimizer of I over this set belongs to X{uq,2j/G) 

The crack predicted by the second model (i.e. V) has no edges. Indeed, suppose 
that V is a smooth manifold. Then the "cut" procedure eliminates the edge if it 
belongs to the sound region. It follows that the edge of the crack lies in a damaged 
region. In other words, around the edge it is at least a region containing a fine 
mixture of damaged and undamaged material. 



I want to thank to Prof. Luigi Ambrosio for his suggestions, especially concerning 
the SBD approach to the first model. 
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